Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 1 February 2008 (MN M?eX style file vl.4) 



A maximum- likelihood approach to removing radio sources 
from SZ observations, with application to Abell 611 

William F. Grainger, Rhiju Das^, Keith Grainge, Michael E. Jones, 
Riidiger Kneissl, G. G. Pooley, Richard Saunders 

Astrophysics, Cavendish Laboratory, Madingley Road, Cambridge CB3 OHE, UK. 



o 
o 

(N 
oo 



> 
oo 

O 



1 February 2008 



ABSTRACT 

We describe a maximum-likelihood technique for the removal of contaminating radio 
sources from interferometric observations of the Sunyaev-Zel'dovich (SZ) effect. This 
technique, based on a simultaneous fit for the radio sources and extended SZ emission, 
is also compared to techniques previously applied to Ryle Telescope observations and 
is found to be robust. The technique is then applied to new observations of the cluster 
Abell 611, and a decrement of — 540± 125 /iJy beam~^ is found. This is combined with 
a ROSAT HRI image and a published ASCA temperature to give an Hq estimate of 
■ Mpc-^ 



52^24 km s-i 
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1 INTRODUCTION 

This paper is concerned with the subtraction of radio sources 
which would otherwise contaminate o r obliterate detections 
of the Sunyaev-Zel'dovich (SZ) effe ct ( ^unyaev fc Zel'dovicl 
|l97[ ; sunyaev & Zel'dovich 1972 ) towards galaxy clusters 



. The work described here is in connection with the Ryle Tele- 
, scope (RT, see e.g. Jones et al. 2001 and Grainge et al. 
» I ■ 2001b), but the issues are relevant to all cm- wavelength SZ 
' observations with interferometers (see e.g. Reese et al. 2000). 
For a massive cluster at moderate or high redshift, the flux 
that the RT detects from the SZ effect at 15 GHz is typically 
—500 /.iJy on its shortest baselines. This is sufficiently faint 
that radio sources will almost invariably be present with 
comparable or greater amplitudes. Thus removing the effects 
of radio sources is an essential step. We describe and com- 
pare two past methods of measuring SZ decrements in the 
presence of sources as well as a maximum-likelihood method. 
We then apply this to new RT observations of cluster A611 
which we combine with X-ray data to estimate Hq. All co- 
ordinates are J2000 and, except where otherwise stated, we 
use an Einstein-de-Sitter world model. 



2 REMOVING RADIO SOURCES FROM SZ 
OBSERVATIONS WITH THE RT 

As the RT is an interferometer with a wide range of base- 
lines, it can simultaneously measure the extended SZ flux 
and the fluxes and positions of the small angular size radio 
sources. Figure |^ illustrates the variation of SZ flux with 
baseline for the RT when observing a massive Abell cluster 
at z = 0.171. The variation with redshift is slight over the 
range z = 0.15-5, as shown in Grainge etal. ( |2001b| ). The SZ 
effect has effectively been completely resolved out for base- 
lines above ~ 1.5 kA, and so these "long" baselines can be 
used to measure sources and so remove their effects from the 
SZ signal seen on the "short" baselines. The measurements 
are simultaneous and, of course, at the same frequency, and 
so the spectral index is unimportant. Variability is unim- 
portant if the t elesco pe conflguration does not change; see 
Grainge et al. (199(3) for details. By choosing an interfer- 
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ometer configuration such that there are more long base- 
lines than short, it is possible to optimise the observations 
to achieve good signal-to-noise for the SZ effect without it 
being dominated by noise from unsubtracted sources. 

There are three methods of source subtraction that have 
been applied to RT data: 

(i) using Clean, which has been used for the bulk of pub- 
lished SZ measurements from the RT (e.g. Grainge et al. 
1996); 

(ii) the matrix method, which was first used by Grainge 
mm-, and 
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Figure 1. The flux density as a funetion of baseline for a splierical 
King-model cluster with central electron density = 10'*m~'^,/3 = 
0.65, core radius = 60", = 7.8 X lO'^ K, z = 0.171 and ^0=50 
km Mpc~^, as observed with the RT. (Prom Grainge 1996). 
The shortest RT physical baselines are 870 A; this becomes some- 
what shorter with projection. 



Table 1. The source fluxes and positions in the simulated data 
set. The convention used in this and other tables is: positive ARA 
is an increase in the RA value, so the source is to the East (to 
the left on conventional maps) and positive ADec is an increase 
in Dec, i.e. to the North. 



Source number 


Flux 


Offset from pointing centre 




(/^Jy beam~^) 


(arcseconds) 


1 


2960 


-10, 10 


2 


910 


35, 15 


3 


255 


-50, -40 


4 


170 


-120, 100 


5 


100 


0, 


6 


80 


60, -60 



Table 2. Fluxes and positions as measured by Maxfit after 
CLEANing with four Clean boxes. The noise on each flux is 38 /iJy. 

Source number Flux Offset from pointing centre 



(/xjy beam ^) (arcseconds) 

1 3040 -9.6, 9.6 

2 970 34.5, 15.8 

3 270 -49.9, -45.8 

4 170 -120.1, 91.9 



(iii) the FluxFitter method, first used by Das (1999) 
and used here in modified form. 



This sequence began with Clean, the classical radio- 
astronomy image deconvolution technique (see e.g. Greisen 
1994 and Perley et al. 1989). The matrix method is a linear 
method for removing the effects of sidelobes, and FluxFit- 
ter - using maximum likelihood - addresses the problem of 
simultaneously fitting both radio sources and SZ decrement. 



2.1 The test data 

The three methods can be explained and compared with an 
example of a simulated dataset containing both sources and 
an SZ effect. The simulation is of a 54 x 12- hour long ob- 
servation of a field at declination 44°. The uv coverage is 
based on a standard Cb configuration for the RT; in this 
configuration four aerials are parked on an east-west rail- 
track at locations 36, 72, 90 and 108 m from the closest 
fixed aerial (see Grainge et al. 1996 for details). The point- 
source fiuxes and positions are shown in Table ^, and the 
noise level was set to 7 mjy visibility"^'^^, corresponding to 
200 /iJy day~^^^, as expected for a standard RT observa- 
tion with five aerials. The SZ decrement is based on Abell 
2218, and is modelled as an isother mal ellipsoid with a King 
electron density profile (King 1972) at the centre of the map 
with a central electron density of 10*m~'^, f3 = 0.65, a tem- 
perature of 8.7 X 10^ K and core radii of 60" and 40" on 
the sky and 49" (= V60 x 40) along the line of sight. The 
central temperature decrement for this cluster is 0.82 mK 
and, as observed with the RT, the cluster gives —660 pJy 
beam~^ on the shortest baseline. 



2.2 The CLEAN method 

For this method, a dirty (i.e. unCLEANed) map of all 
t he baselines longer than 1.5 kA is produced within AlPS 
(Greisen 1994). Natural weighting is used to give the best 



possible signal-to-noise ratio at the expense of larger side- 
lobes. The map is CLEANed in the standard way by placing 
Clean boxes around the obvious sources. After deconvolu- 
tion, the source fluxes are measured, using the AlPS verb 
Maxfit, which interpolates the position and value of the 
maximum fiux density. The measured fiuxes are then re- 
moved from the visibility data using the task UvsuB. 

The sources and the positions found in the simulated 
dataset are listed in Table ^, and the CLEANed map is shown 
in Figure ^. The noise in the map due to the system temper- 
ature is 38 /iJy beam~^. The six model sources are labelled; 
note that only four brightest were found. 

A map made with baselines longer than 1.5 kA was con- 
sistent with noise after these four sources had been sub- 
tracted. The SZ flux observed was —675 ± 72 mJy beam~^, 
at a position (—4, —2) arcsec from the pointing centre. The 
SZ values were measured with Maxfit from a map made 
with the baselines shorter than 1.0 kA. 

A comparison of the sources found (Table ^ and the 
sources actually in the model (Table ^ shows the how well 
the Clean method works. The most glaring problem is that 
only four out of the original six sources have been detected. 
As lower-frequency surveys such as FIRST, NVSS or opti- 
cal images often provide information about the source en- 
virnoment at 15 GHz, the long-baseline map is rarely used 
in isolation. As such, a second test was performed, placing 
the Clean boxes as before, but adding small Clean boxes 
around the locations of the other two weaker sources. The 



Table 3. Fluxes and positions from the Clean method with six 
Clean boxes. 



Source number 

1 
2 

3 
4 

5 



Flux 

(/^Jy beam~^) 

3020 

950 

240 

180 

85 

82 



Offset from pointing centre 

(arcseconds) 

-9.6, 9.5 

34.5, 15.8 

-47.6, -42.0 

-121.6, 98.4 

12,-7 

70,-70 
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2.3 The matrix method 



The matrix method was initally developed by KG (Grainge 



Figure 2. A CLEANed map of the simulated data made with 
basehnes longer than 1.5 kA. Source 1 has had 2500/iJy beam~^ 
of flux already removed to aid identification of weaker sources. 
The contour levels are — 55x(l, ^2, 2) fijj beam~^ (dashed), 
55x(l, ^2, 2, ys, 4, ...) fijy beam~^ (solid). The beamsize 
FWHM (45" X 26" at 4°) is shown in the bottom left. The noise 
in the map duo to the system temperature is 38 ^Jy bearn"^. 



positions and fluxes of the six sources are shown in Table ^ 
The noise on the long-baseline map is 38 fiJy beam~^. Thus 
additional radio or optical information leads to the detection 
of the two faintest sources. 



2. 2. 1 Source finding 

Source finding is done in a non-linear, iterative manner. A 
typical source-subtraction process may involve many itera- 
tions of map making, running the subtraction algorithm on 
the sources found, mapping residuals, finding another source 
and then adding that into the model. 

A high signal-to-noise source is easy to identify, and 
causes no problems. However, the process becomes subjec- 
tive around signal-to-noise ratios of 4 to 3.5. This is the point 
at which the non-Gaussian, correlated statistics in the im- 
age plane conspire with the high RT sidelobes (when using 
a few antennas for source finding) to make source identifi- 
cation more difficult. 

Radio sources in the clusters we observe have angular 
sizes smaller than the maximum resolution used. However, 
the number of Clean components is generally much larger 
than the number of sources in the field. With the test data, 
only four sources were identified, but AlPS produced a model 
with 28 Clean components. Since the real sources are ev- 
idently point sources, this is over-modelling the data; and 
potentially biasing. There is also a degree of subjectivity in 
the placing and sizing of Clean boxes. 



1996), and we here describe it and assess its performance. 
In the matrix method, sources are first identified from a 
unCLEANed long-baseline map, and Maxfit used to mea- 
sure the positions and the fluxes of the sources. The con- 
volution that occurs when observing with an interferometer 
means that the flux on the dirty map of the j"^ source, 
S'dirtyj, is given by 



(1) 



where S'sky,i is the true flux on the sky of the i"* point source, 
Bi^j is a factor due to the synthesised beam that depends on 
the displacement (on the sky) between the source in ques- 
tion and the z"" source, and Pi is the primary (envelope) 
beam attenuation. The value of Bi^j is directly measured 
from the dirty beam produced by HORUS. There is thus a 
matrix equation linking the measured dirty fluxes with the 
true sky fluxes, which is solved by inverting the matrix. This 
method has the advantage over Clean in that sources that 
are measured in the map to be point sources are modelled 
with one flux and position. With the simulated test data, the 
four sources identified with the Clean method were used. 
The resulting matrix was 



/ 2875 \ 




( 1 


0.13 


-0.03 


0.16 \ 


/ Ssky,! 


473 


- 


-0.13 


1 


0.02 


-0.04 1 


'S'sky,2 


105 




-0.03 


0.02 


1 


°f j 


Ssky.S 


V 585 / 




^ 0.16 


-0.04 


0.12 




\ 5'sky,4 



where the vector on the left side of the equation is a mea- 
surement of the dirty flux (in ^Jy beam~^) at each point. 
The value 5sky,i is the beam-attenuated flux on the sky of 
source x, using the same labelling as for sources found in 
the Clean method. 

The matrix method is linear, an apparent advantage 
over the Clean method. If after the first subtraction at- 
tempt some sources are still present, then the additional 
terms for the matrix can then be measured and the solution 
recalculated. 

Table ^ shows the fluxes found. A comparison with the 
sources known to be in the model (Table ^ shows that the 
flux of source three is significantly different. There are two 
reasons for this. Firstly, not all the flux has been subtracted 
- see Figure ^ Secondly, the positions used have been deter- 
mined from the unCLEANed map, by searching for extrema 
with Maxfit. Occasionally, the relative positions of sources 
are such that one source is in the steepest part of a sidelobe 
of another, and Maxfit does not flnd an extremum that 
corresponds to that source. This has occurred in this case 
with source 4. In this situation, Maxfit is not used, and 
the position from the Clean method is used. The flux value 
from the unCLEANed map in that position is then used. The 
reliance on the Clean method combined positional prob- 
lems and with difficulties automating the process for large 
numbers of sources means that the "matrix method" has 
only been applied in simple cases with a few well separated 
sources. 
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Table 5. Source positions and fluxes as reported by FluxFitter 
using the Clean model as an initial guess. The formal error on 
each flux is 30 ^Jy beam~^. 



8^ 43'"" 4l 



20" 



Figure 3. Long - i.e. greater than 1.5 kA - baseline map for the 
test data after the sources in Table ji] have been subtracted. The 
greyscale range is from -200 (light) to 200 fj,3y beam~^. 



Table 4. Source fluxes and positions from the matrix method. 



Source number Flux 

(/ijy beam^^) 

1 2970 

2 877 

3 8 

4 142 



Offset from pointing centre 

(arcseconds) 

-9, 9.6 

34, 15 

-47, -42.8 

-122, 99.9 



2.4 The FLUXFITTER method 

In an attempt to overcome the problems of the Clean 
and matrix methods, the FluxFitter algorithm was in- 
troduced. The algorithm that it uses is straightforward; 

(i) An initial model of the sky is made, using a set 
of parameters that represents the positions and fluxes of 
each source, including the SZ decrement. This model is de- 
termined from the long-baseline RT maps (either raw or 
CLEANed) and, for the SZ parameters, from the X-ray im- 
age. 

(ii) The flux that the RT would observe is calculated for 
every visibility point. 

(iii) The misfit between the model and real uv data is 
then calculated as x^- 

(iv) The parameters are then varied to minimise ■ 

The best-fitting parameter values are then used to sub- 
tract the radio sources; this is done within AlPS. There are 
two advantages to this method. It works almost entirely in 
the aperture plane; only the source identification and ap- 
proximate position finding is done in the image plane. Work- 
ing in the aperture plane is preferred because the noise dis- 
tribution is known to be Gaussian. The second advantage of 
this method is in point (i): a simultaneous fit to the posi- 
tions and fluxes of the point sources and the SZ decrement 
is a clear improvement over either of the previous methods 



Source number 

1 
2 

3 
4 



Flux 

{fj.3j beam"'^) 

3031 

966 

208 

184 



Offset from pointing centre 

(arcseconds) 

-9.5 ±0.1, 9.6 ±0.2 

34.5 ±0.3, 15.3 ± 0.5 

-47.5 ± 1.2, -42.8 ± 2.8 

-122.3 ± 1.6, 98.2 ± 2.8 



as there is no arbitrary "long" and "short" baseline split, 
and it allows full use of all visibility data, which increases 
the signal-to-noise ratio for point-source measurements. 

The initial parameters for the sources are still estimated 
by iterating the Clean method. In a complex situation with 
both bright and faint sources, the map is CLEANed, and then 
positions measured. These sources are then subtracted and 
the subtracted data then mapped again. This loop can be 
performed many times to estimate the number of sources 
and their approximate positions and fluxes. In a less complex 
situation. Clean is not used, and sources are approximately 
subtracted and then the data remapped. Again, this is just 
to provide an initial guess for FluxFitter. 

Currently only the amplitude of the SZ decrement is 
varied. The other parameters that describe the decrement - 
position, core radius and /9-parameter for the cluster - are 
all fixed in advance from the X-ray data. This is done as the 
radio data do not constrain well the core-radius or /3 value; 
with present telescopes, the X-ray measurements constrain 
the core radius, /3 value and position much better. 

FluxFitter was run twice on the simulated test data. 
The fiuxes and positions from the Clean method were used 
as an initial guess for the first run. The fiuxes and posi- 
tions from this run are shown in Table |^ the errors are 
those reported by FluxFitter - see below. After subtrac- 
tion with UvsuB, the map of baselines greater than 1.5 kA 
was consistent with noise. For a second run, the fluxes and 
positions from the Clean method were used, and the po- 
sitions of the additional two sources were also used. The 
fluxes and positions from this run are shown in Table ^ 
This table shows that the two additional sources are de- 
tected with good significance. That the overall noise level of 
30 /iJy beam~^ is lower than that for the Clean method is 
not surprising as all the baselines are being used in the deter- 
mination of the fluxes and positions. After subtracting the 
six reported sources, a map of baselines shorter than 1.0 kA 
shows a decrement of flux —700 ± 65 mjy beam~^ at an 
offset (9", 4"). FluxFitter itself finds a central decrement 
of 0.87 ± 0.08 mK, close to the set value of 0.82 mK. 

As an additional check, a third run of FluxFitter 
was performed, and the four sources found from the Clean 
method and two random points were used as the initial 
guess. In this case, FluxFitter reported that both of the 
two random "sources" had fiuxes below the noise level and 
very large positional-errors. 

FluxFitter also reports error bounds. As Figure ^ 
shows, the contours are elliptical and orientated along 
the variable axes, which shows that the parameters are in- 
dependent. The error on each parameter is calculated by 
finding the parameter values at which the reduced i^" 
creases by 1. 
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Table 6. Source positions and fluxes as reported by FluxFitter 
using additional information. The formal error on each flux is 
30 /ijy beam~^. 

Source number Flux Offset from pointing centre 



(/xjy beam ^) (arcseconds) 

1 3024 -9.9 ±0.1, 9.9 ±0.2 

2 970 34.9 ±0.5, 15.6 ±0.7 

3 212 -47.8 ±1.8, -40.8 ±3.6 

4 155 -121 ±2.6, 98.8 ±4.1 

5 140 7.9 ±3.8, -3.1 ±4.2 

6 132 67.4 ±3.3, -64.8 ±4.6 




-200 -150 -100 -00 00 100 100 200 



Source flUK [liUy trom best tit) 

Figure 4. contours for the flux and position in Dec for the 
brightest source in the test field. The spacing between contour 
levels in each figure is such that the reduced value increases 
by 1 between each contour. 




-10 

Source tlu^ [^Jy trom best tit) 



Figure 5. contours for position for the brightest source in 
the test field. The spacing between contour levels in each figure 
is such that the reduced value increases by 1 between each 
contour. 

The error-bound reporting was checked by simulating a 
point source with differing signal-to-noise ratio. 500 obser- 
vations of a single point source were simulated; the signal- 
to-noise was kept constant for groups of 10 simulations, and 
the position was held constant for all the simulations. The 
visibilities were simulated with Gaussian noise. The known 
position was then fed to FluxFitter as an initial guess, 
and the best-fitting position and flux recorded. It was found 
that the quoted error bar does enclose the position for 67% 
of the simulations. It was also found that the uncertainty of 
the position, that is the size of a, varies as the inverse of 
the signal-to-noise. This result is shown in Figure Note 



o 

o 




1 10 100 

S/N 

Figure 6. Signal-to-noisc ratio verses the uncertainty (in arcsec- 
onds) of the position. Fitting a straight line to this data gives 
uncertainity /arcsec = (13.7 ± 0.5)(signal/noise)~'''^^^''''^'^ 



that this relation holds down to very low signal-to-noise ra- 
tios. The result is useful for determining whether a tentative 
source found with the RT at low signal-to-noise has a po- 
sition coincident with higher significance data, for example 
from NVSS or POSS. 

2.4-1 Possible improvements 

Source recognition is clearly the biggest problem that still 
remains. It is the only step that is still performed in the 
map plane rather than aperture plane. There are computa- 
tional issues involved here: producing a map and identifying 
sources "by hand" is possible and fairly cheap in computer 
time, but the noise in the map plane is non-Gaussian. Min- 
imising the misfit between the data and a given number of 
sources is also cheap; but allowing the number of sources to 
vary vastly increases the complexity of the problem and the 
time required. It is possible that more advanced minimising 
techniques such as simulated annealing (see e.g. Press et al. 
1993) or using massive inference techniques will make this 
possible and robust in the aperture plane. 

2.5 Comparison of results 

Tables ^, ^ and ^ show the results of three different meth- 
ods for source fitting. For both the Clean and FluxFitter 
methods, the results with six sources are considered. The po- 
sitions and fiuxes put into the model are shown in Table |l| 
The resultant position and depth of the SZ decrements af- 
ter subtraction are shown in Table ^. The parameters are 
all measured from dirty maps, made with baselines shorter 
than 1 kA. Note that the matrix method has a less deep SZ 
fiux density, and that the FluxFitter and Clean method 
values are statistically consistent with the expected value for 
this model cluster (—660 /iJy beam~^). Also, the positions 
of the SZ decrement are fully consistent with the model as 
the beam size is around 180". It is not surprising that the 
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Table 7. Parameters for the resulting SZ decrements from the 
simulated data. 



08' 



Subtraction 

method 

Clean 

Matrix 

FluxFitter 



SZ flux 
(/ijy beam~^) 
-675 ± 72 
-500 ±90 
-700 ±65 



position 
(arcsoconds) 
-4, -2 
18, 
9, 4 



flux of the SZ decrement from the matrix method is less 
deep as the central 100 /iJy source has not been subtracted. 
Also, the matrix method does an incomplete subtraction of 
the sources it does find, resulting in more contamination of 
the SZ signal. Note that this does not imply that the matrix 
method will always will give a lower SZ decrement if the 
source subtraction is incomplete. 

All three methods benefit from prior knowledge of the 
source distribution on the sky. This can be estimated from 
looking at lower frequency surveys such as NVSS or FIRST. 
Most falling-spectrum sources, i.e. with a > (where 5* cx 
v~" where v is the observing frequency and a is the spectral 
index) will be detected i n NV SS and or FIRST. Howe ver, a s 
shown by Cooray et al. ( 1998| ) in clusters and Taylor ( 200C ) 
generally, there are rising spectrum sources, i.e. with a < 
that are present at 15 GHz and not detected in NVSS and 
FIRST. 



3 OBSERVATIONS OF ABELL 611 



3.1 X-ray 



AbeU 611 is a cluster ai z ^ 0.288 ( [Crawford et al. 199q ) 
originally identified by Abell ([l957[). It has an 0.1- 2.4 keV 
luminosity of 8.63 x lO""* W ( [Bohringer et al. 2000| ), with a 



temperature of 7.95^0.11 x 10' K ( [White 20001 ). White de 



rived this value from a 57-ks ASCA exposure by considering 
both a single-phase and two-phase cooling model. The tem- 
perature values found for the bulk of the gas are statistically 
equivalent and a mass deposit rate of Olj!,^^ Mq yr~^ was 
found for the cooling model. 

The 17-ks ROSAT HRI observation from April 1996 
is shown in Figure ^ The image contains two bright pix- 
els, which, on comparison with the POSS image, are co- 
incident with a large galaxy. These pixels were removed. 
We calculate an X-ray emissivity constant of 1.29 xlO~^^ 
counts s~^ from 1 m'^ of 7.95 x lO^-K gas of electron den- 
sity 1 m~^ at a luminosity distance of 1 Mpc, assuming a 
metallicity of 0.21 ± 0.07 solar and an absorbing H column 
of 4.88 X 10^" m-^ 

The best-fitting model parameters were /3 — 0.59, core 
radii of 26" and 24" with a position angle of the major axis 
of 101° and a central electron density of no = 11.6 x 10^ 
m~"^ (assuming a core radius along the line of sight of 25 — 
(24 X 26)^''^ arcsec and Ho= 50 km s"^ Mpc"^). There is 
a degeneracy between the core radii fitted and /3 but this 
has no significant effect on Ho (see Grainge et al. 2001b and 
Jones et al. 2001). 



3.2 RT observations 

A611 was observed for 16 sets of 12 hours between Novem- 
ber 1994 and January 1995 with the RT in configuration Gb. 



36" 00 




8 01 20 



8 01 00 



40" 



Figure 7. ROSAT HRI image of A611. The exposure time is 17 
ks. The greyscale range is to 32 counts. 



Table 8. The radio sources in the FIRST catalogue within 400" 
of the RT pointing centre for A611. 

Flux (mJy bcam~^) 
RA Dec Peak 

8 1 20.248 36 5 9.3 1.25 ±0.13 
8 54.948 36 9 6.1 1.10 ±0.13 

Flux and phase calibration and overa ll data reduction strat- 
egy are described in Grainge et al. ( 2001a ). Three days of 
data - taken in bad weather - were rejected after examining 
the 1-day maps and noise levels. A map of the combined 13 
days of data using baselines longer than 1.5 kA had a noise 
level of 70 /iJy beam~^ , and only one source was visible, with 
fiux 299 AtJy beam"^ at RA 8*' 0™ 57=1 Dec. +36° 3' 40". 
This source was removed with UvsuB, using the flux and 
position from the dirty map. A long-baseline map of the 
subtracted data was consistent with noise, with no other 
sources in the fleld. 

Table |^ lists the sources found in the FIRST catalogue 
around the pointing centre for A611. The NVSS catalogue 
contains no sources in this region. Neither of the two FIRST 
sources is detected at 15 GHz and the source that is present 
at 15 GHz is not detected at lower frequencies. 

FluxFitter was then run using the X-ray data to pro- 
vide a model of the SZ decrement and using all the baselines. 
Again, the initial guess was deflned by the 1.5-kA-only flt- 
ting. The source was found to be at RA S*" 0"" 57^ ±0.9 
Dec. + 36° 3' 35"0 ±8 with a flux of 188 ± 65 /xJy beam'S 
which is lower than the long-baseline only values. As the 
angular size of A611 is small, it is likely that the SZ signal 
was contaminating the "long" -baseline map. Note that this 
source is not detected in the FIRST survey, and so in this 
case prior knowledge from a lower frequency survey has not 
helped. Figure ^ shows a GLEANed map of baselines shorter 
than 1 kA after this source has been subtracted. The decre- 
ment (as measured from the map) is — 540±125 /xjy beam"'^ 
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Figure 8. The SZ effect in A611. The contour levels are -480, 
-360, -240 and -120 /ijy beam"! (solid) and 120 /ijy beam"! 
(dashed). The map has been CLEANed; the restoring beam, which 
is 92" by 350" FWHM at a position angle 3.4°, is shown in the 
bottom left. 



at RA 8^ 0™ 57"3 Dec. + 36° 2' 38". This location is 3" 
in RA and 36" in Dec away from the X-ray cluster location. 
Considering the Clean beam used is 92" x 350", this is a 
good positional agreement between the X-ray and SZ obser- 
vations. The slight extension to the south is not significant. 



3.3 Ho determination 

From the source subtracted dataset and the X-ray parame- 
ters of the cluster, it is po ssible to estimate Ho, as described 
in Grainge et al. ( 2001h| ). The likelihoods for different Ho 
values are calculated and are shown in Figure 0. The best- 
fit Ho for A611 is then 52tll km s"^ Mpc"^ The error 
quoted is due to the noise in the SZ measurement, and does 
not include any of the other sources of error in the determi- 
nation. The addit ional s ources of error are described fully 
in Grainge et al. (2001b). For A611, the error from the SZ 
measurement is by far the most important, and the final Ho 
value is 52tjg km s~^ Mpc~^ for an Einstein-de-Sitter world 
model. Assuming a world model with = 0.7, Qm = 0.3, 
Ho=59t.i^ km s~^ Mpc~^ from this cluster. 



4 CONCLUSIONS 

The problem of radio source contamination in interferomet- 
ric SZ observations and methods to remove it have been 
investigated, demonstrating the following. 

(1) The non-linear Clean method can work well, but 
does not use all the available information and can over- 
complicate the problem. 

(2) The matrix method, though linear, fails in typical sit- 
uations such as the one simulated here. The failure is mainly 



20 40 60 80 100 

Hubble Constant (knn/s/Mpc) 

Figure 9. Likelihood plot for different Ho values from fitting to 
the source-subtracted SZ data from A611. 



due to the high sidelobes from the Ryle Telescope; they make 
it difficult to determine accurate positions, and so fluxes, for 
sources close to each other on a map. 

(3) FluxFitter uses all the available information and 
produces the simplest model. It solves simultaneously for 
sources and SZ decrement, and it works with with the vis- 
ibilities, where noise is known to be Gaussian, rather than 
in the map plane where the noise is correlated. 

(4) All three techniques suffer from the problem of source 
identification, which is currently performed in the image 
plane where the noise characteristics are complex. Source 
identification can be aided by prior information, for exam- 
ple from lower-frequency surveys. 

(5) The positional uncertainty, as determined in the 
aperture plane, is found to vary as uncertainity/arcsec cx 
(signal/noise)~°'^^^° °^ even at signal-to-noise ratios below 
nominal detection limits. The constant of proportionality 
will be a function of the interferometer used. 

(6) Observations of the cluster A611 with the Ryle Tele- 
scope give a 4.3-cr detection of an SZ decrement, and combi- 
nation with X-ray data gives and estimate of Ho~ 521^ jg km 
s~^ Mpc~^, assuming an Einstein-de-Sitter cosmology, and 



59 



+27 



km s ^ Mpc ^ using = 0.7 and fl-^ 



0.3. 
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